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We investigate the large distance behavior of the electric and magnetic propagators of hot SU(2) 
gauge theory in different gauges using lattice simulations of the full 4d theory and the effective, di- 
mensionally reduced 3d theory. Comparison of the 3d and 4d propagators suggests that dimensional 
reduction works surprisingly well down to the temperature T = 2Tc- Within statistical uncertainty 
the electric screening mass is found to be gauge independent. The magnetic propagator, on the 
other hand, exhibits a complicated gauge dependent structure at low momentum. 



I. INTRODUCTION 



Understanding the large distance behavior of the static electric and magnetic propagators in hot SU(N) gauge 
theory is of interest for several reasons. First of all it is related to the screening phenomenon in a hot non-Abelian 
plasma. In fact, the concepts of electric and magnetic masses extracted from these propagators form to a large extent 
the basis for our intuitive understanding of screening in non-abelian gauge theories Furthermore, the concept of 
screened electric propagators finds application in refined perturbative calculations (hard thermal loop resummation) 
. Non-perturbative calculations of these propagators thus can provide a bridge between perturbative and non- 
perturbative descriptions of the electric screening phenomenon. 

Of course, a draw back of such an approach is that the gluon propagator itself is a gauge dependent quantity. One 
thus has to question to what extent the results extracted from these propagators have a physical meaning. The poles 
of the gluon propagator at finite temperature were proven to be gauge invariant in perturbation theory However, 
for static quantities like the Debye mass static magnetic fields can lead to a breakdown of perturbation theory. The 
Debye mass cannot be defined perturbatively beyond leading order 1^,0. Non-perturbatively the problem of gauge 
independence of the screening masses thus is an open question. 

It has been shown that at high temperature the large distance behavior of a SU(N) gauge theory can be described 
in terms of the dimensionally reduced effective theory, i.e. the three dimensional adjoint Higgs model |P-p^. The 
screening masses extracted from gauge invariant correlators were studied in terms of the effective three dimensional 
(3d) theory p"2[-Jl7[| and were compared with the corresponding results from four dimensional (4d) simulations p8| , p^ . 
The relation between the propagator masses and the masses extracted from gauge invariant correlators was discussed 
in 1^ . The screening masses extracted from gauge invariant correlators correspond to masses of some bound states of 
the 3d effective theory and are several times larger than masses extracted from propagators. 

In the present paper we extend our earlier studies on the electric and magnetic propagators . Contrary to 

Refs. [p|, pOpl| ] where propagators were studied only in Landau gauge we consider here a class of generalized Landau 
gauges, the Coulomb gauge, the maximally Abelian gauge and the static time averaged Landau gauge. We study the 
propagators in terms of the 3d effective theory as well as in the 4d theory at finite temperature. A comparison of the 
propagators in the full 4d and in the 3d effective theory provides further evidence for the applicability of dimensional 
reduction. A detailed study of finite size effects in the propagators shows that the picture of magnetic screening given 
in Refs. |^,^|T] needs to be revised to some extent. First results from our analysis of the magnetic sector have been 
presented in Ref. p^ . 

The paper is organized as follows. In section II we discuss questions related to gauge fixing and the definition 
of propagators in finite temperature SU(2) theory as well as in the dimensionally reduced effective theory, the 3d 
adjoint Higgs model. Section III contains our main results on electric and magnetic propagators and the analysis 
of their gauge and volume dependence. Finally we give our conclusion in section IV. In an Appendix details of the 
determination of the parameters of the effective theory are discussed. 
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II. GLUON PROPAGATORS IN FINITE TEMPERATURE SU(2) GAUGE THEORY 



In this section we define the actions we use for our simulations in three and four dimensions and introduce our basic 
notation for gauge fields, gluon propagators and the different gauges we have analyzed. 



A. Actions in three and four dimensions 



In four dimensions (4d) all our calculations are performed with the standard Wilson action for SU(2) lattice gauge 
theory, 



(1) 



where U^{x) £ SU{2) are the usual link variables and f34 — ^/gl- In three dimensions (3d) the standard dimensional 
reduction process leads us to consider the 3d adjoint Higgs model 



i + \h\ \t,AI{-s.)+x(\t,AI{-s.) 



(2) 



where 13^ now is related to the dimensionful 3d gauge coupling and the lattice spacing a, i.e. (3^ = 4 /ffq a. The adjoint 
Higgs field is parameterized by hermitian matrices Ao = """^o the usual Pauli matrices) iQ. Furthermore, 

X parameterizes the quartic self coupling of the Higgs field and h denotes the bare Higgs mass squared. The relation 
between the 3d and the 4d couplings will be discussed in an Appendix. We also note that the indices /i, v, of course, 
run from to 3 in four dimensions and from 1 to 3 in three dimensions. Although we will not always mention this 
difference explicitly it should be obvious from the context how various sums that appear have to be interpreted. 



B. Gauge fixing 

As we want to analyze properties of the gluon propagator, which is a gauge dependent quantity, we have to fix a 
gauge on each configuration on which we want to calculate this observable. In the past most studies of the gluon 
propagator have been performed in Landau gauge. Here we will consider a class of A-gauges which are generalizations 
of gauges that have been introduced in to smoothly interpolate between Landau and Coulomb gauge. In the 

continuum these gauges correspond to the gauge condition 

^A^a^A^^O . (3) 

On the lattice these gauges are realized by maximizing the quantity 

TrEA^[/^(a;) . (4) 

In 4d the Landau gauge condition thus corresponds to the case A^ = 1 for all ^ — 0, ..,3, while the Coulomb gauge is 
given by Aq = and A^ = 1 for i = 1, 2, 3. In the latter case we have to impose an additional gauge condition for 
the residual gauge degree of freedom. We do this by demanding 

E C/o(a;o,x) = Wo , (5) 

X 

to be independent of xq. In addition to the A-gauges we also consider the Maximally Abelian gauge (MAG) which 
can be realized by maximizing the quantity p^ ] 
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^Tr[a3[/^(x)a3(7t(a;)] , (6) 

with (T3 being a Pauli matrix. Also in this case one has to fix a residual gauge degree of freedom which we do by 
imposing a U(l)-Landau gauge condition JE6[. In 4d SU(2) gauge theory we also consider the static time-averaged 
Landau gauge (STALG) introduced in p7| , p8| [. In the continuum it is defined by 

3 

9oAo(a;o,x) = 0, Y.Y.^-^^ = (^) 

On lattice this gauge is realized in two steps. First we maximize the quantity Tr^^ Uq{x). In the second step we 
maximize the quantity 

i—l X 

performing xq -independent gauge transformations. 

While the notion of Landau and Maximally Abelian gauges carries over easily to the 3d case we have to explain a 
bit more in detail our notion of A-gauges in 3d. We have considered two versions of A-gauge: 

Ag-gauge: diAi + d2A2 + XsdaAa = (9) 



Ai-gauge : XidiAi + ^2^2 + dsA^ = (10) 



The Ai-gauges are more closely related to the 4d A-gauges considered by us; in both cases the rotational symmetry 
of the gauge condition is broken in a direction orthogonal to the X3 (or z) direction which we are going to use for 
separating the sources in our correlation functions. Furthermore, we introduce in 3d the so-called Coulomb gauges 
(ci gauges) which fix the gauge in a plane transverse to the z-direction, 

ci -gauge : cidiAi + ^2^2 = . (11) 



Of course, as in the 4d case we again need an additional gauge condition for the residual gauge degree of freedom. 
The case ci = 1.0 is the usual Coulomb gauge. This gauge has the advantage that a positive definite transfer matrix 
exists in the z-direction. 

In our numerical calculations the gauge fixing was performed using a standard iterative algorithm accelerated by 



an overrelaxation step 1 29 1 as well as by using a stochastic overrelaxation algorithm |30| . 

When analyzing observables in a fixed gauge one also has to address the question to what extent Gribov copies 
can influence the result. This problem will not be discussed here. It previously has been studied in 4d SU(2) gauge 
theories at zero temperature |3l[ as well as at non-zero temperature [p0| . In both cases no influence of Gribov copies 
on the gluon propagators was found within the statistical accuracy achieved in these studies. 

C. Gauge fields and gluon propagators 

Definitions of the lattice gauge field ^^(x) utilize the naive continuum relation between the lattice link variables 
U^{x) and the gauge field variables in the continuum limit, i.e. Uf^{x) — exp (igaAf^{x)). A straightforward definition 
thus is 

M^)^7^ [U^.{x) - Ulix)] , (12) 
2iga 
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which in the continuum Hmit differs from the continuum gauge fields by 0{{ag)^) corrections. Here a is the lattice 
spacing and g either the 3d or 4d bare gauge coupling. Other possible definitions for A^(a;) which formally lead to 
smaller discretization errors were considered in [ ^2| . It was found there that up to an overall multiplicative constant 
the difference in magnetic and electric propagators due to different definitions of the gauge field is much smaller than 
the statistical errors. We note that the definition ( p^ of the gauge field assumes that the link variable C/^ is close to 
the unit matrix. This seems to be the case for all gauges considered except the Coulomb gauge. In the case of Coulomb 
gauge the gauge fixing procedure does not force the temporal link to be close to the unit matrix. Only the spatial links 
turn out to be close to the unit matrix. For instance, for links averaged over a lattice volume of size 12^ x 24 x 4 we find 
for /34 = 2.512 (corresponding to 2Tc) < TtUq/2 0.576(4) and < >= 0.90089(6). These numbers should 

be compared with the corresponding ones in Landau gauge: < TiUo/2 >= 0.8879(1) and < Trt7i/2 >= 0.88709(7). 
One can expect that TiUo/2 gets closer to 1 as one approaches the continuum limit. Indeed for the Coulomb gauge 
on a 24^ X 48 X 8 lattice at /3 = 2.74 (also corresponding to 2Tc) we find < Tr;7o/2 >= 0.625(6). However, from these 
numbers it is clear that very large values of Nt, i.e. large couplings /34, are necessary to get a meaningful definition 
of the Aq field. Therefore we will use the 4d Coulomb gauge only to analyze magnetic propagators. 
The lattice gluon propagators in d-dimcnsion are defined as 

\ X3,b I 

Here Q^^ixz) is a sum over all gauge fields in a hyperplane orthogonal to X3, 

Qt(^3) = a'^-^E4(^^'^3), (14) 

with J^^Jyx) = ^Ti[Af^{x)a^]; x± ~ {xi,X2) in 3d and x± = {xo,xi,X2) in 4d, respectively. Furthermore, fl = N'^Nz 
is the lattice volume in 3d and il = N'^NzNt in 4d 0. The electric and magnetic propagators are then defined in the 
usual way (see e.g. Refs. pO| , ^ ), 

DE{z)^Da^{z) , Dm{z) = \{D^i{z) + D22{z)) . (15) 

£'33(2) is not included in the definition of Dm{z) because it is constant in A and Coulomb gauges. We will also 
consider momentum space propagators which are obtained through a one dimensional Fourier transformation, 

^MM(fc) = « ^''"'^A.^W'^ith fc= -J^, n = 0,l,....7V,/2 (16) 



We use the standard definition for the momentum space magnetic propagator (see e.g. |33| , |34| ) 

3 

3 



1 ^ 

DMik) = -Y,D^^ik). (17) 

We note that we include here D33 in the definition of the momentum space propagators in order to take into account 
the contribution of the constant component 1)33(2), which only infiucnccs the zero mode contribution D^^{k = 0). 
The electric propagator in momentum space is simply defined by DE{k) = Doo(fc). In order to absorb additional cut- 
off effects in momentum space propagators we find it usefull to analyze these in terms of the momenta p = |2 sin(fc/2)| 
rather than the lattice momenta k. In the following we thus will use p instead of fc as a definition for our lattice 
momenta. For the analysis of the long distance behavior of these correlation functions it is customary and, in fact, 
quite instructive to consider local masses, which are defined by 

A(^) cosh(m,(z)(z-iVj2)) 

(iSj 



D,{z + 1) cosh(m,(2)(z + l-7V,/2))' 



^We will also use the notation V — N'^Nz for the spatial volume in 3d and 4d. 
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with i = E, M. If the propagators decay exponentially starting from some value of z, the corresponding local masses 
will reach a plateau. 

All our 4d simulations have been performed at a temperature T — 2Tc. The values for the 4d coupling cor- 



responding to this temperature were taken from Ref. 
lattice volumes used in our 4d simulations. 



n 



and are given in Table I together with the corresponding 



/34 = 2.512 


8^ X 16 X 4 
12^ X 24 X 4 
16^ X 32 X 4 
20^ X 40 X 4 
24^ X 48 X 4 
28^ X 56 X 4 
322 X 64 X 4 


/3i = 2.740 


16^ X 32 X 8 
242 X 48 X 8 



TABLE I. Couplings and lattice volumes used in the 4d simulations of finite temperature SU(2) gauge theory. 



We give the choice of parameters for our 3d simulations in the next subsection. The cutoff dependence is discussed 
for the magnetic and electric propagators separately in section III. 



D. 3d versus 4d calculations 



Most of the results we are going to discuss in the following section have been obtained through simulations in the 
dimensionally reduced version of the 4d SU(2) gauge theory at finite temperature. As we do want to compare our 
results obtained in 3d with corresponding results in the 4d theory obtained at a temperature T = 2Tc we should check 
that dimensional reduction yields reliable results at temperatures this close to the critical point. In Refs. ll^-|l^] 
it was shown that the effective theory is capable to describe the long distance behavior of some gauge invariant 
correlators in SU(2) and SU(3) gauge theories at these temperatures. A similar conclusion has been reached for the 
2d dimensionally reduced version of the 3d SU(3) gauge theory at finite temperature 

We have performed a detailed analysis of propagators calculated in Landau gauge in 4d and 3d at T = 2Tc. The 
latter has been simulated at values of h corresponding to the metastable region of the 3d adjoint Higgs model, i.e. for 
f] — 11^ h — —0.3773 and x = 0.099 (see appendix for details). The 4d simulations have been performed a,t (3 — 2.74 
on lattices given in Table I. In order to compare the propagators obtained from simulations in 3d and 4d we have 
normalized them to unity at distance zT = 0. The results are shown in Figure |l|. As one can see from the figure the 
3d and 4d results for the electric as well as for the magnetic propagators are in excellent agreement. It is interesting 
to note that good agreement between 3d and 4d propagators is found already at relatively short distances, although 
the 3d effective theory is expected to describe the 4d physics only at distance zT > 1. The above result implies that 
dimensional reduction works quite well even at T = 2Tc and can be established by comparing gauge fixed observables. 
In Figure || we compare the electric and magnetic propagators calculated from full 4d and effective 3d theories at 
2Tc for the maximally Abelian gauge. As one can see from the figure a good agreement between 4d and 3d results 
exists also here, although it is a priori not clear in which sense the 3d Maximally Abelian gauge corresponds to the 4d 
Maximally Abelian gauge. The agreement between the magnetic propagators in Figure |^ shows that the two gauges 
are quite similar. A possible reason for this agreement is the decoupling of the Aq field from the dynamics of spatial 
(magnetic) gauge fields which will be discussed in the next section. We note, however, already here that the magnetic 
propagators calculated in Landau and Maximally Abelian gauges, respectively, show quite a different long distance 
behavior. This is evident from the comparison of Figure 1 and Figure 2 and will be analyzed in much more detail in 
the next section. Also the apparent volume dependence visible in these figures will be discussed later in more detail. 
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FIG. 1. Comparison of 4d and 3d data for the electric (left) and magnetic (right) gluon propagator in Landau gauge calculated 
at T = 2Tc. 
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FIG. 2. Comparison of 4d and 3d data for the electric (left) and magnetic (right) gluon propagator in Maximally Abelian 
gauge calculated at T = 2Tc. 



III. NUMERICAL RESULTS ON ELECTRIC AND MAGNETIC GLUON PROPAGATORS 

In this section we will discuss our numerical results obtained from simulations in three and four dimensions. The 
main purpose of this investigation was to quantify the gauge dependence of the electric and magnetic gluon propagators 
and analyze to what extent gauge invariant masses can be extracted from the long distance behavior of the gluon 
correlation functions. 

In Figure ^ we show results from a calculation of the magnetic propagator in 4d finite temperature gauge theory 
in different A-gauges including the Coulomb gauge limit. All the A-gauges yield identical magnetic propagators. In 
fact, one should expect that all A-gauges do lead to identical propagators on an infinite spatial lattice. For static 
configurations {doAo = 0) the gauge condition for the A-gauges, XSqAq + diAi = is up to multiplicative factors 
identical to the 3d Landau gauge condition. As we have shown that the 4d propagators in Landau gauge can be 
mapped onto the 3d propagators this has to be the case also for the 4d propagators calculated in A-gauges. The same 
holds for STALG. 

We also note, that magnetic propagators show strong volume dependence at distances zT>\. This effect may be 
traced back to the influence of zero mode contributions to the propagators which are different in different gauges 
and give a volume dependent positive contribution to the correlation functions. Zero mode contributions are most 
prominent in the long distance behavior of the correlation functions calculated on finite lattices. 
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FIG. 3. The magnetic propagator in A gauges on a 16^ x 32 x 8 lattice (left), as well as on 24^ x 48 x 8 and 12^ x 24 x 4 



lattices (right) at T = 2Tc. Here also the propagator in STALG is shown. The propagators were normalized to 1 at 



0. 



In Figure ^ we summarize our results for electric propagators. As in the case of magnetic propagators, electric 
propagators calculated in different 4d A-gauges as well as in STALG agree with each other as expected. Contrary to 
magnetic propagators the electric propagators show no significant volume dependence. We also note that the cut-off 
effects seem to be small both for electric and magnetic propagators (cf. Figures |and|). The issues of volume and 
cut-off dependence will be discussed separately for electric and magnetic propagators in the following two subsections. 
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FIG. 4. The electric propagators in different gauges. Shown is the electric propagator measured on 16^ X 32 x 8 (left) and 
on 24^ X 48 X 8 and 12^ x 24 x 4 (right). The propagators were normalized to 1 at z = 0. 



A. The magnetic gluon propagator 



In this section we will discuss in more detail our results on magnetic propagators. Previous lattice calculations of 
the magnetic propagators in hot SU(2) gauge theory in 4d ||2|] and 3d |g] gave evidence for its exponential decay in 
coordinate space and thus indicated the existence of a magnetic mass. A non-zero magnetic mass was also found in 
analytical approaches based on gap equation - 39 1 and a non-perturbative analysis of 2 -I- 1 dimensional gauge theory 
pO| . Nonetheless, it has been questioned whether a magnetic mass in non-abelian gauge theories does exist [41 4^. In 
Ref. the Landau gauge propagators of 3d pure gauge theory were studied on fairly large lattices and were found 
to be infrared suppressed 0. Such a behavior clearly is in conflict with the existence of a simple pole mass. In order 



^Similar results were found in 4d SU(N) gauge theory at zero temperature p^ , ^ . 
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to further investigate this problem simulations on larger lattices are needed to explore the long distance regime of the 
correlation functions. This can be achieved in the 3d reduced theory which gives a good description of the 4d theory 
even at temperatures a few times T^. We thus concentrate, in the following on a discussion of results obtained from 
our calculations in 3d. In fact, it also has been observed already in earlier calculations that the magnetic propagators 
of the 3d adjoint Higgs model are very close to the corresponding propagators of 3d pure gauge theory. We thus 
further restrict our analysis of the magnetic propagators to the limit of 3d pure gauge theory. Where appropriate we 
will perform a comparison with results obtained from the 3d adjoint Higgs model and the 4d SU(2) gauge theory. 

In order to get control over the propagators also in the continuum as well as in the infinite volume limit we have 
performed calculations for = 5, 5.5, 6, 8 and 16 on different lattice volumes. The simulation parameters are 
summarized in Table II. 
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48^ X 64 




563 








643 
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TABLE II: Lattice volumes used at different values of in our 3d calculations of the magnetic gluon propagator. 



It has been pointed out in |45[ | that zero modes can give sizeable contributions to gluon correlation functions. In 
terms of the lattice gauge fields the zero mode contribution is defined as the expectation value of the average gauge 
field, 

'^m = ^E^m(^) ■ (19) 

X 

The zero mode contribution is apparent in the momentum space propagators where one finds from cq. |l^ for vanishing 
momentum 

D,,{p = Q)^a''nJ2{{r,f) ■ (20) 

a 

By taking the inverse Fourier transformation of Eg. ([l^) one can easily see that the zero modes give a constant 
contribution {Nza)~^ Df^p_{jp ~ 0) to the coordinate space propagator D^^{z). In fact, such a contribution does 
qualitatively explain the overall volume dependence of the magnetic propagators observed in our calculation. The 
contribution from zero momentum mode fluctuations is expected to vanish in the infinite volume limit. We thus 
expect that at fixed distance zT the propagators approach their asymptotic, infinite volume values from above as the 
lattice size is increased. This is indeed the case, as can be seen in Figure ^ where we show the 4d and 3d magnetic 
propagators calculated in Landau and Maximally Abelian gauges on different size lattices at T = 2Tc. As one can see 
from the figure the volume dependence is quite different in these two gauges. In the case of Landau gauge we observe 
that the magnetic propagators in coordinate space decay faster with increasing lattice size. In the Maximally Abelian 
gauge, however, the propagator does not exhibit any sizeable finite size dependence. 
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FIG. 5. Volume dependence of the magnetic propagator in coordinate space calculated in Landau (left) and Maximally 
Abelian (right) gauges. The propagators were normalized to 1 at 2; = 0. Shown are results from 3d simulations at /^a = 8 
and 4d simulations, which are compared at similar values of the physical volume in units of . These volumes correspond to 
lattices of size 16^ x 32 x 8, 24^ x 48 x 8 and 16^ x 32 x 4 in 4d and to 16^ x 32, 24^ x 48, 32^ x 64 and 48^ x 64 in 3d. The Pa 
values for our 4d simulations are given in Table I. 
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FIG. 6. Large distance behaviour of the magnetic propagator in coordinate space calculated in Landau gauge. The propa- 
gators were normalized to 1 at z = 0. 

For volumes yr'^>300 the magnetic propagator calculated in Landau gauge becomes negative for zT>2. In Figure 
^ we show the large distance behavior of the coordinate space propagator calculated in 4d and 3d on lattices with 
similar 3d lattice volume. A similar behavior of the coordinate space propagator was found in other A-gauges. The 
propagators calculated in MAG, however, stay positive for all lattice sizes and distances accessible in our calculation, 
i.e. up to zT ^ 5. 

As was mentioned above the magnetic propagators in the 3d adjoint Higgs model are not very sensitive to the pres- 
ence of the adjoint Higgs field. This is illustrated in Figure |^ where we compare magnetic propagators in momentum 
space calculated in the 3d adjoint Higgs model with corresponding results obtained in the pure gauge theory at the 
same value of f]^- This also shows that the effect of the adjoint Higgs field on the magnetic propagator decreases with 
increasing temperature. In fact, for temperature T>10Tj, we find no visible effect of the adjoint Higgs field on the 
magnetic propagators. In what follows we will therefore mainly discuss the magnetic propagators in the limit of pure 
gauge theory. 
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FIG. 7. Momentum space magnetic propagator obtained from 3d pure gauge theory and the adjoint Higgs model with 
couplings corresponding to temperatures 2Tc and 3Tc. Shown are the magnetic propagators in pure gauge theory at /Js = 8 (a), 
in 3d adjoint Higgs model ai /Ss = 8, x — 0.09, h = —0.4846 corresponding to the temperature 3Tc (b), in pure gauge theory 
a,t 13 — 5.5 (c) and in 3d adjoint Higgs model at /3 = 5.5, x — 0.099, h = —0.7528 corresponding to the 2Tc (d). 

In Figure |^ we show the momentum space magnetic propagators in Landau gauge obtained from 3d pure gauge 
theory at /Ss = 5 and Ps — 8. For momenta p > O.65I the propagators in units of g| are volume and /33-independent. 
For small momenta, however, they are strongly volume and also dependent. One can clearly see that for large 
volumes the propagators reach a maximum at non-zero momentum and start decreasing with decreasing momenta. 
This is a direct consequence of the negative propagators in coordinate space found in Landau gauge for zT>2. 

The volume dependence of the magnetic propagator in momentum space is strongest at p = 0. On the other hand at 
Ps = 5 one can see that the value of the propagator in the vicinity of its maximum is essentially volume independent 
for the three largest volumes. A similar behavior of the momentum space propagator was found in other A-gauges. 
This is demonstrated by Figure Note that although the propagator is a gauge dependent quantity the position of 
the peak seems to be gauge independent in the class of A-gauges considered here. 

The existing rigorous bound on the infrared behavior of the gluon propagator in momentum space imply that it is 
less singular than p~^ in 4d and p^^ in 3d It has been further argued that it is likely to vanish at zero momentum 
in the thermodynamic limit |47|. T hese considerations were also extended to the case of magnetic propagators in gauge 
theories at finite temperature |[48[ . We have analyzed this question by performing studies of the volume dependence 
of the zero momentum propagator. The data for Dm{0) are shown in Figure |l^ for different P3 values. For = 5 
and ^ 8 the data have been fitted to the ansatz 

DMiO)gl = aiVgl)-^+b . (21) 

We have performed fits with b — and b > 0. In the first case we have obtained rather large values (typically the 
value of x^/d.o.f was between 4 and 10) for both P^ = 5 and P3 = 8. The three parameter fits give a reasonably 
good for both P3 values. Using the whole range of volumes one finds a = 51(11), & = 2.06(15) and z = 0.34(3) 
with x^/d-o.f. = 1.8 for /Jg = 5 and a = 52(9), = 2.08(33) and z = 0.31(3) with with x^/d-o.f. = 1.8 for P3 = 8. 
Although the volume dependence of Dm{0) seems to be /Jg-independent, the value of the exponent z as well as the 
value of Z?A/(0) in the infinite volume limit strongly depends on the range of volumes used in the fit. For example 
fitting the data on Dm{0) for P3 = 5 using lattice volumes from 30'^ to 96^ one gets : a — 281(139), — 2.52(56) 
and z = 0.55(5). In order to firmly determine the functional form of the volume dependence of -Dj\/(0) and the value 
^m(O) in the infinite volume limit simulation on even larger lattices are necessary. 

In order to confirm the infrared suppression of Dm{p) and the existence of a maximum in Dm{p) at non-zero 
momentum we have analyzed in more detail the interplay between infinite volume and continuum limit. 
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FIG. 8. The momentum space magnetic propagator at /Js = 5 (open symbols) and — 8 (filled symbols) calculated in 
Landau gauge. 
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FIG. 9. The momentum space magnetic propagator at /Ja = 5 on 72 lattice in different A-gauges. 





P=5 . 






P=8 ^ 
=16 ^ 


— 1 
— 8 — 1 






$ /: 

M - 








n 


, "M 



































1e-06 1e-05 0.0001 0.001 

1/(Vg3^) 



FIG. 10. Volume dependence of Dm{0) at different values of fig. The lines represent fits with the ansatz given in Eq.(pi| 
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We have seen that the propagators are infrared suppressed for large enough volumes at /33 = 5 and 8. However, the 
magnitude of the propagator at small momenta shows a clear /33-dependence. Therefore, it is important to estimate the 
behavior of the magnetic propagator in the continuum limit. We have seen that a.t f]^ — 5 the magnetic propagator 
on a 60'^ lattice is already close enough to the corresponding infinite volume limit. Therefore we have performed 
additional calculations of the magnetic propagators at (3^ — 6 and on 72'^. The lattice volumes 60"^, 72^ and 96"^ at 
(33 — 5, 6 and 8 correspond to the same physical volume. Moreover, for small momenta p < 0.5g|, values of the 
propagators at approximately same values of momenta are available. In this region of momenta we have performed an 
extrapolation to the continuum limit fitting data at difi'erent /Js with the ansatz a + b/ (is. The result of this analysis 
is summarized in Figure 1 1 , 
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FIG. 11. The Landau magnetic propagator at different values of (3^ and the continuum extrapolation. 

As one can see from this figure the general structure of the propagator, in particular the existence of a maximum 
at non-zero momentum is preserved in the continuum limit, although the volume dependence of the momentum space 
propagator rapidly increases with decreasing momentum for p/g^ < 0.5. 



B. The electric gluon propagator 



Let us now turn to a discussion of the electric gluon propagator. The volume dependence of the electric propagator 
turns out to be quite different from the magnetic ones. The electric propagators show exponential decay at large 
distances in all gauges considered. 

In Figure ^ the Landau gauge electric propagators in momentum space are shown for T — 2Tc and T — 9200T'c- In 
the latter case we performed the analysis only in the 3d effective theory at x = 0.03, /3 = 16 and h — —0.2085. As one 
can see from the figure the volume dependence of the electric propagators is indeed quite different from the volume 
dependence of the magnetic propagators. In particular the electric propagators show no sign of infrared suppression. 
For fixed lattice geometry, i.e. fixed ratio N/Nz the volume dependence of the zero mode contribution seems to 
approach as expected for the propagator of a massive particle. Moreover, the infinite volume limit is almost 
reached on the largest lattices used in our simulations. We find a similar volume dependence also in other gauges. 

Since we are interested in extracting the screening mass from the electric propagator it is also important to address 
the question of volume dependence of the electric mass. In Figure |l^ the local electric masses in Landau gauge are 
shown for different lattice volumes for (3^ — 16, h — —0.2085 and x — 0.03. As one can see from the figure, the value 
of the plateau (which in fact determines the screening mass) is essentially volume independent, although there is some 
volume dependence in the local masses at short distances. Thus the screening masses can be estimated even from the 
smallest lattice for these values of the parameters. We will use this fact in the discussion of the gauge dependence of 
the electric screening mass. 
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FIG. 12. The Landau gauge electric propagator in momentum space at 2Tc (left) and T = 9200rc (right). 
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FIG. 13. Local electric masses in Landau gauge for /3 — 16, h = —0.2085 and x = 0.03. The solid line is the value of the 
electric mass obtained from an exponential fit of the electric propagator calculated on a 32^ x 96 lattice. The dashed lines 
indicate the uncertainty in its value. 

Since the electric propagators show exponential decay at large distances in all gauges considered it is natural to 
ask whether the electric masses extracted from them are gauge independent. There is a formal proof that poles of 
finite temperature propagators are gauge independent to any order of perturbation theory ||^. However, as discussed 
in section I the Debye mass is not calculable in perturbation theory beyond leading order. Therefore it is not clear 
whether the arguments of Ref. |^ apply for electric screening mass. In terms of the effective 3d theory the screening 
of static electric fields is related to propagation of the adjoint scalar field. Since the 3d theory is confining the adjoint 
scalar field is not a physical state ( the physical states are some bound states) and therefore there is no physical 
principle which guarantees its mass to be gauge invariant. Guided by the analogous problem in 3d scalar QED it 
was conjectured in Ref. that the large distance behavior of the high temperature non-Abelian electric propagator 
is dominated by a brunch cut singularity. Although this is different from the simple pole obtained in leading order 
perturbation theory this still leads to an exponential decay of the electric propagators in coordinate space with a gauge 
invariant screening mass. Thus it is clear that also the question of gauge (in)dependence of the electric screening mass 
is non-trivial in general. 

We have studied the gauge dependence of the electric propagators for two sets of parameters of the 3d effective 
theory x = 0.03 , h = —0.2085, = 16 and x = 0.03, h = —0.1510, = 24 both corresponding to the temperature 
T = 9200rc- In this investigation lattices with spatial volume 16^ x 32, 24^ x 48 and 32^ x 96 were used. From studies 
of the Landau gauge propagators we know that electric masses for /?3 = 16 can be reliably extracted already from 
16^ X 32 lattice. We have investigated electric propagators in different gauges introduced in section II. B. For = 16 
most simulation were performed on 16^ x 32 lattice. For calculations with A3 gauges we also used a 32^ x 96 lattice. 
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For = 24 propagators in all gauges except the Maximally Abelian gauge were measured on 32^ x 96 lattice. Due to 
the large number of iterations required for fixing the Maximally Abelian gauge propagators were measured only on a 
242 X 48 lattice at P3 = 24. 




FIG. 14. The local electric masses calculated on 32^ x 96 for x = 0.03, h = -0.1510 and 03 = 24 in A3- gauges, in 3d 
Coulomb gauge (ci — 1) and in Maximally Abelian gauge. The solid line is the value of the electric mass obtained from an 
exponential fit of the electric propagator in Landau gauge calculated on a 32^ X 96 lattice for f3 = 16, h ^ -0.2085 and x = 0.03. 
The dashed lines indicate the uncertainty in its value. 

Our main results for — 24 are summarized in Figure |lj, where local electric screening masses in Coulomb, 
Maximally Abelian and A3- gauges are shown. As one can see from the figure local masses are strongly gauge 
dependent at short distances. For zT>2 they, however approach a plateau which is gauge independent within the 
statistical accuracy reached in our calculation. 

We also have performed simulations on a 16^ x 32 lattice for (3^ = 16 using Coulomb gauges with ci — 0.1, 10, 
as well as Ai -gauges with Ai = 0.1, 10 and Maximally Abelian gauge. The corresponding results for local electric 
masses are shown in Figure |l5|. Again within statistical accuracy no gauge dependence of the plateau of the local 
electric masses can be observed, though the local masses have quite different z-dependence. 
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FIG. 15. The local electric masses calculated for x — 0.03, h — —0.2085 and P3 = 16 on 16^^ x 32 lattice. Shown are the 
local electric masses calculated in ci gauges for ci =0.1 (a) and ci = 10 (b), in Ai gauges for Ai =0.1 (c) and Ai = 10 (d) and 
Maximally Abelian gauge (e). The solid line is the value of the electric mass obtained from an exponential fit of the electric 
propagator in Landau gauge calculated on a 32^ x 96 lattice (3 = 16, h = -0.2085 and x = 0.03. The dashed lines indicate the 
uncertainty in its value. 
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IV. CONCLUSIONS 



We have investigated electric and magnetic gluon propagators in the high temperature phase of SU(2) pure gauge 
theory in different gauges. The propagators were calculated directly in full 4d theory and also in 3d effective theory. 
It was shown that the effective theory can describe the electric and magnetic propagators remarkably well down 
to temperature T = 2Tc. We find that electric propagator can safely be extrapolated to the infinite volume limit, 
within the statistical accuracy of our present analysis and the class of gauges considered here its long distance 
behaviour is gauge independent and yields a gauge independent electric screening mass which is compatible with 
earlier determinations. 

The magnetic propagator, however, has a complicated volume and gauge dependent infrared structure which is 
not compatible with a simple pole mass. While it still shows a simple exponential decay at large distances in the 
Maximally Abelian gauge it starts getting negative for zT>2 in a class of A-gauges which includes the Landau gauge. 
This leads to the infrared suppression of the propagator in momentum space and is incompatible with a simple pole 
in the magnetic gluon propagator as it was deduced from earlier studies [p| j2C|pl| , |36[ | which were limited to shorter 
distances and smaller lattices. Nonetheless, we find in all gauges that magnetic correlators are screened at large 
distances. 
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APPENDIX 



We will discuss here the procedure we followed to fix the couplings of the 3d adjoint SU(2) Higgs models. Two 
different approaches were proposed for fixing the parameters appearing in (|^), the usual perturbative dimensional 
reduction |^,^ and non-perturbative matching analyzed in Here we have used both approaches. 

The lattice gauge coupling f]^ is related to the dimensionful 3d continuum gauge coupling by the standard 
relation 

* = ^- P^) 

This relation basically determines the lattice spacing a in terms of the dimensionful coupling g|. The Higgs self- 
coupling and the 3d gauge couplings are related to the renormalized gauge coupling g{p-) of SU(2) gauge theory in 
MS scheme. The corresponding relations were calculated in Ref. ||l^ to 2-loop level. 

To fix the temperature scale one needs to choose the renormalization scale /i which fixes the temperature dependence 
of the 4d gauge coupling. For this we will use information on the temperature dependence of the spatial string tension. 
The analysis of Ref. |5^] showed that the spatial string tension can be well described in terms of 3d effective theory. 
Furthermore it was observed that the temperature dependence of the spatial string tension of finite temperature 
SU(2) gauge theory can well be described by a simple formula 0.334(14)(7^(r)T with g{T) given by the 1-loop 
renormalization group relation. The string tension of the pure 3d SU(2) gauge theory was found to be 0.335(2)g^ [|2). 
The string tension of the 3d adjoint SU(2) theory has been calculated in jlj]. Within statistical errors it turned out 
to be independent of the scalar couplings and equal to the string tension of the pure SU(2) gauge theory. Based on 
these observations we have fitted the data on the spatial string tension from using 0.335{2)g^ with 1-loop level 
couphng and the relation = 1.06A-jg^ from |21]. The fit yields ^ = 18.86T. Wc have also fitted the data for the 
spatial string tension using the 2-loop formula for g^ from p^ ]. However, it turned out that the simple 1-loop formula 
fits the 4d data on spatial string tension much better. Based on these observation the 3d gauge coupling gl and the 
Higgs self-coupling were chosen according to the 1-loop version of formulas (2.13) and (2.15) from 

9l^9'{f^)T, (23) 
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with fi = 18.86r. In Ref. Q the parameters and x were chosen according 2-loop formulas from |12| and ^ = 2nT. 
We have checked that differences between the present choice and that used in Q are about 10%. In Ref. (l4| the 
parameters of the effective theory were also calculated using the 2-loop formulas from ||l^, however, the 4d gauge 
coupling constant was determined from 1-loop formula in MS'-scheme with fi ^ 7.0555. In Ref. the screening 
masses extracted from gauge invariant correlators were calculated in the effective 3d theory and compared with the 
corresponding results from 4d simulations ||l8|,|9| . A rather good agreement between the screening masses calculated 
in 3d effective and full 4d theories was found. However, the 3d gauge coupling constants in |Q were considerably 
larger than ours which would overestimate the spatial string tension by 30%. Nevertheless this contradiction can 
easily be resolved by noticing that the set of parameters corresponding to 5Tc in ||l^, namely x = 0.104 and y = 0.242 
would roughly correspond to the temperature 2Tc in our approach. Using the results from for these values of 
X and y and the value of g|(2Tc) — 2.89T we get for the two smallest screening mass the values to(O^) — 2.90(3) 
and m(OI) ~ 3.91(8) (for details about classification of different gauge invariant screening masses see Ref. [Q). 
These numbers should be compared with those obtained in 4d simulation in ||l^ at T = 2Tc, m(O^) = 3.06(12) and 
m(OI) = 4.06(12). Thus it seems that with our present choice of the parameters the effective 3d theory could describe 
all static quantities measured so far in lattice simulations. We also note that in the case of 2-|-l-dimensional gauge 
theory at finite temperature it has been shown that the spatial string tension and gauge invariant screening masses 
can be simultaneously well described by the effective 2d theory In 2+1 dimensions the situation is simplified by 
the fact, that the gauge coupling has no renormalization scale dependence due to superrenormalizability of the theory. 

Generally we have used the non-perturbatively determined values for ft, However, in some cases the propagators 
were also calculated for values of ft given by perturbative dimensional reduction. The bare Higgs mass parameter h 



could be related to the mass parameter y = m'jjQ/g^ of the continuum adjoint Higgs model |12 

h='^y- 3-1759114(4 +5x) _ 1 / ^^^^ _ 3 ^ ^ ^ ^ \ ^^^^ 

with rriDo being the tree- level (from the point of view of the effective theory) Debye mass. At 1-loop level one has: 



y-^-w^- (26) 

From Eqs. ( p5| ) and (^^ we have calculated values of h corresponding to the 1-loop level dimensional reduction. In 
this case simulation were performed in the metastable brunch of the broken phase (see Ref. ||l^ for details). 

Numerical values for the couplings of the 3d effective action used in most of our calculations are summarized in 
Table HI together with the corresponding values of T/Tc- For other values of P3 used in our calculations we have 
calculated the coupling ft from Eq. (E5|) using the y values from Table III. 
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0.2274 
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-0.2611 


0.3556 


-0.2622 


0.4007 


-0.2700 


0.2755 


0.070 
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16 


-0.2528 


0.4408 


-0.2490 


0.5009 


-0.2588 


0.3470 


0.050 


72.1 


16 


-0.2365 


0.5914 


-0.2314 


0.6721 


-0.2437 


0.4756 


0.030 


9200 


16 


-0.2085 


0.9279 


-0.2006 


1.0544 


-0.2181 


0.7758 



TABLE HI . Numerical values for the coupling ft and continuum mass parameter y = rnj^Q/g'^ of the effective 3d 
theory in our simulations and the corresponding values of the temperature. I and II denote two sets of ft-values 
taken from Ref . |q| . The last two columns give the coupling ft and mass parameter y obtained from the perturbative 
relations Eqs. ( |25[ ) and (p6|). 

The electric propagators show exponential decay at large distances. Therefore it is possible to extract screening 
masses from them. The procedure of extracting the screening masses from the propagators was described in Ref. . 
In Figure |l^ we compare the electric screening masses extracted from the Landau gauge propagators in the 3d adjoint 
Higgs model calculated on 32^ x 64 and 32^ x 96 lattices with the corresponding results from 4d SU(2) gauge theory 
pl[. The parameters of the effective theory were chosen according to Table HI. 
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FIG. 16. The electric screening masses in units ol the temperature. Shown are the electric masses niE for the first (filled 
circles) and the second (open circles) set of h. The line represent the fit for the temperature dependence of the electric mass 
from 4d simulations. The open triangles are the values of the electric mass for h values obtained from perturbative reduction 
in the metastable region (the last column in Table III). Some data points at the temperature T ~ 70Tc and T ~ QOOOTc have 
been shifted in the temperature scale for better visibility. 

As one can see from the figure the agreement between the masses obtained from 4d and 3d simulation is rather 
good. The electric mass shows some dependence on h. For relatively low temperatures (T < SOTc) the best agreement 
with the 4d data for the electric mass is obtained for values of h corresponding to 2-loop dimensional reduction and 
lying in the metastable region. This fact motivated our choice of the parameters of the effective theory at T = 2Tc in 
section III. For higher temperatures, however, practically no distinction can be made between the three choices of h. 

Before closing the discussion on the choice of the parameters of the effective 3d theory let us compare our procedure 
of fixing the parameters of the effective theory with that proposed in [^. In Ref. the gauge coupling was fixed 
by matching the data on Polyakov loop correlators determined in lattice simulation to the corresponding value 
calculated in lattice perturabtion theory. The resulting gauge couplings turned out to be considerably smaller than the 
corresponding ones used in our analysis, e.g. for T = 2Tc it gave g"^ = 1.43 while our procedure gives g^{2Tc) — 2.89. 
The scalar couplings were fixed according to 1-Ioop perturbation theory ^ . Using this procedure the authors of Ref. 
1^ obtained a good description of the Polyakov loop correlator in terms of the 3d effective theory, however, the spatial 
Wilson loop whose large distance behavior determines the spatial string tension was not well described in the reduced 
3d theory. The reason for this is the fact that the value of the Polyakov loop correlator is cut-off dependent and 
therefore it is not very useful for extracting the renormalized coupling. 
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